Diffusing non-local inflation: 
Solving the fleld equations as an initial value problem 
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There has been considerable recent interest in solving non-local equations of motion which contain 
an infinite number of derivatives. Here, focusing on infiation, we review how the problem can be 
reformulated as the question of finding solutions to a difi^usion-like partial differential equation with 
non-linear boundary conditions. Moreover, we show that this diffusion-like equation, and hence the 
non-local equations, can be solved as an initial value problem once non-trivial initial data consistent 
with the boundary conditions is found. This is done by considering linearised equations about any 
field value, for which we show that obtaining solutions using the diffusion-like equation is equivalent 
to solving a local but infinite field cosmology. These local fields are shown to consist of at most two 
canonically normalized or phantom fields together with an infinite number of quintoms. We then 
numerically solve the diffusion-like equation for the full non-linear case for two string field theory 
motivated models. 

PACS numbers: 



I. INTRODUCTION 

Non-local field equations which contain infinitely many 
derivatives have recently received considerable attention 
in the context of cosmology Q, 0, i, i, S, i, 0, i, B [l^l ■ 
In particular, the effective actions of a number of mod- 
els derived from string field theories have been studied. 
Examples are the ground state tachyon in bosonic cubic 
string field theory fBCSFT) fill and cubic super-string 
field theory (CSSFT) (see [I4] for reviews), and inp-adic 
string theory [ll| . 

Keeping these particular examples in mind, wc will for 
the majority of this paper keep the model we consider 
as general as possible so that our methodology can be 
applied widely. Moreover, though we study only single 
field models it is clear that our method can easily be ex- 
tended to encompass additional fields. As we will discuss 
later, the method can also be generalized to solve cosmo- 
logical perturbations in non-local models, and this is a 
major motivation for its development. The key ingredi- 
ent we require is that the non-local operator which the 
field equations contain be of the form exp(— aD) where a 
is a constant and □ the d'Alembertian operator. As we 
will show, an operator of this form implies that finding a 
solution to the non-local field equations is equivalent to 
finding a solution to a partial differential equation (pde) 
which resembles the diffusion equation. This idea is not 
new (see for example related work in Refs. [1, [1, [l^). 
and has recently been employed in papers very closely 
connected to this work As we discuss in more detail 
later, this recent work solves the diffusion-like pde as a 
boundary value problem, that is, it specifies the initial 
and final state of the field and seeks solutions which in- 



terpolate between these states. Our interest, however, is 
in how the field equations, and hence the diffusion-like 
pde, can be solved as an initial value problem, specifying 
only initial conditions for the field's evolution. 

For linear models an initial value formulation does not 
require us to consider the diffusion-like equation directly, 
since, as we will review, a correspondence between the 
non-local field equations with one field and local field 
equations with an infinite number of fields can be made. 
Thus the problem of specifying initial data for the solu- 
tion of the non-local field equations is reduced to specify- 
ing initial conditions for these local fields |,i (also see 
[15|). For an isotropic cosmology these initial conditions 
correspond to the initial value and velocity for each field 
(we should note that not all fields need to play an active 
role - some can simply be set to zero). In our work we 
extend the understanding of this correspondence with an 
infinite number of local fields by generalizing the linear 
models previously considered with quadratic terms, to 
include potentials with linear terms as well. Moreover, 
we show that if we require real solutions, the fields can 
be interpreted as containing at most two canonically nor- 
malized fields or phantom fields, together with an infinite 
number of quintoms. Let us review the definition of these 
fields. Here and throughout we use the signature of the 
metric (—,+,+,+). A field cp coupled to gravity with 
the action 
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represents a canonically normalized field when s — 1, 
but a phantom when s — —1. A quintom is defined as a 
combination of one canonically normalized field and one 
phantom field interacting through a shared potential and 
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has the action 



II. NON-LOCAL LAGRANGIANS 



We begin by considering a non-local action for a single 
scalar field "0 minimally coupled to gravity, of the form: 
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Ml 



R 



(2) 

Phantom fields violate the weak energy condition (the 
name together with a major industry began with [l7|), 
and hence their energy density can increase with time. 
Quintoms, on the other hand, have the property that 
their energy density oscillates with time. Pairs like these 
where put forward as a simple way to explain the ap- 
parent transition of the equation of state of the universe 
from p -I- P > at large redshifts to p -I- P < at small 
redshift and were coined as "quintoms" [l^. Phantom 
fields raise concerns because naively they are unstable 
to vacuum decay (see for example [l^ for constraints), 
though the question here is complicated as they enter at 
an effective level. For the purposes of this paper, which 
is to determine how to solve the equations of motion, we 
will leave the issue of physicality of phantoms aside. One 
should also be aware of the possibility of classical instabil- 
ities (see Ref. [1^, though the construction in this work 
requires a finite number of derivatives in the theory). 



The major step forward made in our work is to show 
that, for linear models, solving the non-local equations 
using the correspondence technique is completely equiv- 
alent to solving the equations using the diffusion-like pde. 
Furthermore, we determine what initial data for the pde 
is equivalent to introducing each of the local fields. By 
linearising a non-linear model about a particular field 
value, we will be able to use these results to specify ini- 
tial data for non-linear models, which we then evolve into 
the full non-linear regime using the diffusion-like equa- 
tion. Moreover, the understanding we gain from the lin- 
ear models allows us to give an interpretation of the ini- 
tial conditions for the pde, even in the non-linear case, 
in terms of which fields are present initially. 



S = dt^X\J~g 



Mil . / 1 



where □ represents the d'Alembertian operator 
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and in this work we will take 



F{U) = -{l + A^'^aU)e 



(3) 
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which covers a wide range of theories appearing in the 
literature including the p-adic string, the bosonic string 
field theory (BCSFT) and the cubic superstring field the- 
ory (CSSFT). 7'* and a are constants, fixed by the 
particular theory of interest, later we will see a number 
of concrete examples but for now we can treat them as 
free parameters. We should note that in the string field 
theory context our coupling to gravity has been intro- 
duced essentially by hand (as it has been in all other 
works on non-local cosmology) as the minimal covariant 
generalisation of the flat space Lagrangian. 

The scalar field equation which follows from the action 
® is 



(f -I- 4e'^an)e-"°V 
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while the Einstein field equations take the usual form 



G -^T 



(7) 



with Gfj,i, being the Einstein tensor, and the non-local 
stress energy tensor ( first derived in Ref. [f6,]) is given 
by 



Our paper is structured as follows: In section II we in- 
troduce the non-local model we will be considering. Sec- 
tion III reviews the difficultly in solving non-local field 
equations as an initial value problem before introduc- 
ing the diffusion-like pde. Section IV considers linear 
models and their correspondence with an infinite field 
local cosmology, and section V and VI shows how this 
understanding leads to suitable initial data for solving 
the diffusion-like equation. Numerical solutions are pre- 
sented in section VII for two examples. In section VIII 
we discuss how our method can be extended to solve per- 
turbation equations in non-local cosmology, and finally in 
section IX we conclude and outline future directions. 
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III. SOLVING AS AN INITIAL VALUE 
PROBLEM 



In order to solve a scalar field equation with a large 
but finite number of derivatives, for example n deriva- 
tives of Tp, we would need to specify an initial value for 
-0, V", d^^-^tp/dt"^^, which can be given independently 
of one another. For a system with an infinite number of 
derivatives, however, this would require an infinite num- 
ber of initial conditions and moreover it would imply that 
in principle any solution is possible. This argument fol- 
lows since freedom to specify each of the infinite number 
of initial values is equivalent to the freedom in specify- 
ing a full Taylor series solution for the field's evolution 
(assuming of course that the solution is analytic in the 
neighbourhood of the initial conditions) . This is not the 
case, however, since the scalar field equation itself im- 
poses an infinite number of constraints on the allowed ini- 
tial values of the derivatives (as is discussed at length in 
Ref. [21]). For example, for a model with V{ip) = — V''*/4 
and = 0, let us take some initial value ipo — ifj{t — 0) 
and fix the initial conditions such that all derivatives are 
zero except Dilj\t=o. Equation (O then impHes the con- 
dition — an?/)|o = ipo ~ V'o, however, it can be verified, 
in Minkowski space, that operating twice on both sides 
of Eq. ^ by the d'Alembertian operator implies that 
either □V'lo = or -00 = which is incompatible with 
the specified initial conditions. This simple example il- 
lustrates that differentiating the equation of motion im- 
poses additional constraints upon the initial data. In 
general, repeatedly applying the d'Alembertian operator 
to both sides of the scalar field equation generates an in- 
finite number of constraints on the initial data. Our ini- 
tial data cannot therefore be arbitrary, and attempting 
to specify initial data which satisfies all the constraints 
seems hopeless. This is an important point to which we 
will refer later. 



A. The diffusion-like equation 



It seems that the naive formulation of the initial value 
problem for non-local field equations is intractable. We 
will now show that we can reformulate the problem in 
an alternative way by recalling a connection between the 
non-local field equations and a diffusion-like pde [1, 
1,113. Although our discussion can be easily extended 
to an arbitrary form of the metric, we will at this point 
specialize to the homogeneous and isotropic Friedmann- 
Robertson- Walker metric, as our primary interest is non- 
local cosmology. In this case, is a function only of time, 
and the d'Alembertian operator becomes 



□ 
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where H is the Hubble rate, H = d/a, given by the usual 
Friedmann equation 
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where for our non-local theory, p is given by the stress 
energy tensor ^ such that p — —Too, and takes the form 
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a dr 
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and the Raychaudhuri equation also takes the usual form 
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where the pressure P is given by 
7^A 



P = 
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To make the connection with a diffusion-like pde, we 
introduce an auxiliary variable r, and define a new field 
as 



*(i,r) = e-™°0(t) , 



(14) 



which is a function of both time and the auxiliary vari- 
able r. By differentiating ^>{t,r) with respect to r, it is 
simple to see that ^I^ must satisfy the diffusion-like partial 
differential equation 



1 a*(i,r) 
a dr 



(15) 



Moreover, we can see that *(i,l) = e-"°V(i) and 
^'(^,0) — ip{t), hence, in terms of our new field ^(t,r) 
the scalar field equation ([6]) simply becomes 



^{t, 1) - 



9*(i,r) 



dr 



dV{'^{t,0)) 
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(16) 



It is clear, therefore, that if \E'(r, t) is a solution to the 
pde p5|) which satisfies the boundary condition given by 
Eq. (fTH]) at all times, then ip{t) = *(t,r = 0), is a so- 
lution to the non-local field equations, and moreover, H 
can also be determined through Eq. (|10p (and in turn 



the behavior of the scale factor can be obtained). A 
compUcating factor is that the Hubble rate H appears 
in the d'Alembertian operator □, so i? and ^'(t, r) must 
be solved for simultaneously. This is however straight- 
forward, since using the definition of \E'(r, t) and the dif- 
fusion equation (fT5|) . the stress energy tensor ^ can 
be written purely in terms of '^{t,r) for r in the range 
< r < 1, with no non-local operators appearing explic- 
itly. For the FRW case, the energy density, for example, 
is given purely in terms of ^ and its first time derivative 
as 

p = t^[^,{t,0)-^(t,l) + 2V{-^{t,Q)) 



- 4C^a^'(t,0)«'(t, 1) 



dr 




(17) 
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FIG. 1: In the upper diagram we illustrate the method of 
solving the diffusion-like equation as a boundary value prob- 
lem. In the lower diagram, we show how the diffusion like 
equation is solved as an initial value problem. 



B. Solving the diffusion-like pde as a boundary 
value problem 

One method of numerically solving the diffusion-like 
equation which gives rise to a solution to the non-local 
field equations has been explored in previous studies Q. 
In these works the pde (fT5|) is solved as a boundary value 
problem. This approach has the draw back of being more 
restrictive in the solutions it can explore than an initial 
value formulation, but nevertheless it has been highly 
successful and it is well worth reviewing before we pro- 
ceed to the initial value formulation. 

The method attempts to find a solution to the pde 
in the region < r < 1 and —oo < t < oo, visual- 
ized in the upper panel of Fig. [TJ The method requires 
'ii{t — — oo,r), and ^'(i ~ oo,r) to be specified (i.e. the 
boundaries to be fixed). The only obvious consistent data 
for these boundaries is for 5" to be constant for all values 
of r and equal to field values corresponding to turning 
points of the full potential part of the action ([3]), that 
is V + ip'^ /2. At these points all the infinite derivatives 
can be consistently set to zero, and the field is station- 
ary at the maximum or minimum of its potential. The 
other two boundaries ^(t,r = 0) and '^{t,r — 1) are 
then assumed to be related through the condition given 
by Eq. p6p . Then a solution is found using some form 
of relaxation/iterative procedure. The simplest example 
would be as follows: initial trial functions for ^'(t, 0) and 
H{t) for — CX3 < i < CX3 are set, making sure that both 
asymptote to constants at i = — cxd and t = oo which 
are consistent with the field being at an extremum of its 
potential asymptotically in the past and future. The ini- 
tial form of ^'(t, 0) is then used as an initial condition 
for numerical integration of the pde in the r direction 
from r = to r = 1. Once the initial trial function is 



integrated to r = 1, a new initial function is determined 
by employing the boundary condition pB]) , the new trial 
function given by 



*(t,0)r 



/ dV 



Va«'(i,o) 



4? 



dr 



'J'(t,l)oid,(18) 



where {dV/d'i'{t,0))~^{x) denotes the inverse function 
of dV/d'if{t,0). A new trial function for H{t) is also 
determined by employing ^'(i,r)oid in Eq- or in 

Eq. (flO)) . This procedure is then repeated until conver- 
gence is achieved. Obviously variations on this procedure 
can be used, and more sophisticated relaxation or other 
techniques could be employed (also see ^22,] and references 
therein for finding solutions iteratively not directly using 
the diffusion equation in the Minkowski space setting). 



C. Solving the diffusion equation as an initial value 
problem 

Despite its success, the boundary value formulation 
can only produce solutions with the particular initial and 
final conditions where all derivatives of the field tend to 
zero, and hence can only find solutions which evolve from 
one extremum of the field's potential to another. Indeed, 
although interesting and physically relevant solutions can 
be obtained using this method, these requirements are 
rather restrictive. For example, we may wish to consider 
what happens to a non-local field which is initially at 
the maximum of its potential and then subjected to a 
small fiuctuation and ask: what is the full range of pos- 
sible evolutions that can result? Moreover, as we discuss 
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in section IVIIIl the boundary formulation when apphed 
to perturbative equations is not able to produce physi- 
cally interesting solutions. For a perturbation we would 
normally like to specify some initial conditions, for a lo- 
cal field this is usually that the field initially be in the 
Bunch-Davis vacuum far inside the cosmological horizon, 
integrate the system forward in time and determine how 
the perturbation evolves. This is not possible within the 
framework of the boundary value formulation. 

An interesting question, therefore, is: can we solve a 
non-local cosmology as an initial value problem using the 
diffusion-like pde? Such a possibility would allow us to 
specify initial conditions beyond the restrictive ones dis- 
cussed above, and help us to understand the freedom 
allowed in specifying initial conditions and hence the va- 
riety in the possible solutions to the field equations. This 
is the purpose of the present work. 

The idea is illustrated in the lower panel of Fig. [1] and 
we now describe it in detail. We specify some initial 
data {^(^i, r), 4/(^1, r)}, which in turn gives us the ini- 
tial value H{ti) through Eq. (|10p . We then use these 
initial functions together with the differential equation 
(|f Sp and the boundary condition (|16[) in the r direction 
to evolve \['(ii,r) one time step forwards in time. The 
new r) is then used to calculate a new H{t), and the 
procedure is repeated. Continuing this procedure indefi- 
nitely ought to allow us to evolve the initial data forwards 
in time. Naively everything seems as it should be, i.e., 
we specify the two initial conditions we require at the ini- 
tial time boundary (the correct number as the equation 
is second order in time), and have the required condition 
relating the r = and r = 1 boundary. This boundary 
condition is unusual, but numerically its implementation 
is not a problem. There are however two serious issues 
with this idea. The first is that the system is ill-posed 
for generic initial data. This means that high frequency 
modes present in the initial function grow faster than 
lower frequency modes, and this leads to a numerical in- 
stability as any small numerical inaccuracies will show 
up as the introduction of high frequency modes which 
rapidly grow and swamp the solution we are looking for. 
This can be seen clearly for the case of constant H where 
the diffusion equation can be solved exactly for linearised 
boundary conditions. The most general solution is given 
in Eq. (j54p . and here we simply note that by studying 
that solution it can readily be seen that as the solution 
becomes more oscillatory in the r direction, so its growth 
in time increases. It can also be seen that the Hubble rate 
H has a damping effect decreasing the growth rate, but 
that sufhciently oscillatory solutions will always have a 
growing mode. This seems a terminal problem, but we 
will see that in many circumstances of interest the insta- 
bility can be suppressed by regularization methods and 
reliable solutions can still be found. 

The second issue is that the initial data \E'(ii,r) and 
^(ti,r) cannot be arbitrary, but must also satisfy the 
boundary condition p6p . This is a more stringent con- 
dition than one might initially think as this condition 



implies an infinite number of other conditions on the rela- 
tionship between the function ^'(ii, r) evaluated at r = 1 
and at r = 0. These conditions result from successively 
differentiating p^ . Let us then illustrate this for the case 
when ^(ti,r) = which is a consistent initial condition, 
and for convenience we take f ^ = and ^^("0) — — ?A^/4. 
Applying □ successively to the boundary condition (IK 
we find 



dr 



Qj,2 



3^'(ii,0)2 
3^'(ti,0)2 
+18*(ti,0) 



9*(ii,r) 



dr 



Qj.2 



dr 



r=0 
2 



r=0 



(19) 



for the first three operations, where we have used the 
pde itself to replace □ with derivatives with respect to r. 
We can of course continue this process indefinitely. We 
now need to find a function of r which satisfies an infinite 
set of constraints on its derivatives at its boundaries. In 
fact, this infinite set of conditions is simply a reflection 
of the infinite number of constraints that would have to 
be imposed on the initial conditions for -0, -0, etc, which 
we discussed in section Hill The task of finding this func- 
tion 5'(ti,r) is difficult, and we can only hope to find a 
function that can approximately satisfy the conditions. 
We propose to do this by computing the linear evolution 
of -0 about a chosen value A such that 0(i) = A-\- (j){t) 
and hence ^{t,r) = exp(-ran)(A -|- 0(t)) = A + <f>(t, r). 
We will see that this expansion will allow us to find the 
profile in r of ^(ti,r), the problem being much simpler 
because the equation of motion is linearised. This profile 
in turn gives us a profile for 4'(ti,r) approximately. Let 
us then study the linear non-local equations first. 



IV. LINEAR AND LINEARISED EQUATIONS 

In this section we consider the situation when the equa- 
tion of motion ([6]) is linear or it can be linearised by ex- 
panding -0 around a chosen point A such that i(j = A + (f>. 
Our main purpose is to aid our study of non-linear equa- 
tions by linearising them, but it should be noted that 
if the model was intrinsically linear this section furthers 
our understanding of these models in there own right. 
Inserting ip = A + 4> into the action ^ and expanding 
to quadratic order in 0, it results that the equation of 
motion remains unchanged in terms of 0, 



but with a potential ?7(0) given by 



dU{4>) 



U{(j)) = Uo-c(li--m''4>\ 



(20) 



(21) 
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where = V{A) + 



-V'iA) - A, 



-V" (A), a prime means differentiation with respect to ip 
and the quantities are evaiuated at ip = A. Let us now 
define the field <I>(i, r) — e~''"'^0 such that now we have 
the diffusion-hke equation 



1 9$ 

a dr ' 



with the hnear boundary condition 
"9$(i,r) 



dr 



TO^$(i,0) 



(22) 



(23) 



The diffusion-hke equation can be solved by separation 
of variables using 



$(i,r) = 0(t).g(r) + h{r) . 



(24) 



Substituting into Eq. (p2)) we obtain that the functions g 
and h must be of the form g = e"'^'^ h^b e^'^'Vw^ ^ 
while (f) satisfies the local equation 

□</> = -wV-&, (25) 

where k, uP' and 6 are constants. By imposing $(t,0) — 
(/){t) we obtain k — —hjuP' and by making use of the 
boundary condition we obtain h = coJ^ /{m? — 1) 
and the characteristic equation 



'(l-4C^aw2 



(26) 



We note at this point that to solve for (j) it is only re- 
quired to solve a local equation with parameters and 
b. In fact, as it has been done in previous literature, we 
could have postulated l[25|) . substituted it into the equa- 
tion of motion (PO)) and we would have obtained the same 
expression for b and the characteristic equation (|26[) with- 
out any mention of the diffusion-like equation. Indeed, 
the diffusion-like equation is only required when we move 
to the non-linear boundary conditions as in that case we 
cannot localize the field. This section, however, serves 
to illustrate that the two methods are equivalent in the 
linear case. 

Let us consider in detail the roots of the characteristic 
equation ([^5]) . In the case £^ = 0, we have for ^ 0: 



— Ini 



± 2mT i . 



and for 



= 0: 



(27) 



(28) 



where n is an integer. As we will see below, this means 
that in this case the dynamics of a single, massive, non- 
local field can be rewritten in terms of local field equa- 
tions for at most one real field together with an infinite 
set of complex fields. We will see later that the complex 
fields can be coupled with their respective complex con- 
jugates to form two physical real fields which make up a 



"quintom". For a free field, = 0, the dynamics does 
not admit a local description. 

For the case f ^ ^ 0, one obtains for ^ 0: 



where x 



auj — ——r + LambertW (n, x) , 



'exp(-l/4^2)/4^2^ g^^^ fQj. ^2 ^ Q. 



1 

4^ 



(29) 



(30) 



and n is again an integer. Hence, provided 7^ 0, 
a real root exists for x > —0.37 when n — and for 
—0.37 < a; < when n = — 1, otherwise, the roots are 
complex. As we will see, it follows that the dynamics al- 
lows for a description in terms of at most two real fields 
and an infinite number of complex fields. For linearised 
models, the approximate local description resulting from 
an expansion at a given point A in the potential is very 
useful as it gives us an understanding about the possible 
evolution of the field, for example, near an extremum of 
the potential. However, to understand how the solutions 
evolve beyond any linearised regime, or how they inter- 
polate between say a maximum and a minimum one is 
required to go beyond the linear expansion into the non- 
linear regime of the evolution. To do this is one of the 
main aims of our work and we show numerical results for 
two string models in Sec. IVIII 



A. The single real field solution 

Before moving to the non-linear case, it is instructive 
to see how in the linear case a correspondence can be 
made between the non-local dynamics and local dynam- 
ics in terms of canonically normalized fields, phantom 
fields and quintoms, which we defined in the introduc- 
tion. First we will discuss the situation when only a 
single real root to the characteristic equation Eq. P5|) 
is considered. We are interested in real solutions to the 
non-local equations, and so we will only be interested in 
situations in which (jj is real. Inserting ip — A + (j3 into the 
expression for the stress energy tensor ([5]) and employing 
Eq. (|25p , the integrals can now be evaluated and we find 
that eq. ([5]) can be written in terms of a real rescaled 
field 

VP=|B|0, (31) 
with 52 = 74 Aa [1 - 4^2(2 4. iaoj'^)] e""' such that 

TlT' = S9,u Qa^'/^aV + W{^)^ - sd^^d,^ , (32) 

where s is the sign of and the potential for the canon- 
ically normalized field is 



W{ip) = Wo- ship - s^ijj'^Lp^ , 



(33) 
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where b = \B\b, and 



where now the scalar potential involving all the fields is 



Wo = j^X 



1 + auj^ - e""' 



(34) 

Comparing these equations to Eq. ([T|) , we see that when 
is positive ip represents a canonically normalized field, 
while when it is negative it represents a phantom field. 
It is clear, therefore, that is a solution to the non-local 
field equation if 1^9 is a solution to the local Einstein 
and scalar field equations 



Dip 



—uj ip — b , 

^local 

Af2 



(35) 



This remarkable simplification has been explored by a 
number of authors (for example Ref. [1, [j, M , the 
ansatz Eq. ((25|) was first proposed in Ref. [ifl])- In 
particular, in Ref. 6] a general non-local system with a 
hilltop potential of the form (PT|) . but with c = 0, was 
rewritten in the manner presented here in terms of a lo- 
cal field and it was found that in order to act as an in- 
fiationary potential satisfying the current observational 
constraints, the parameters must be Wq « 3.7 x 10~^° 
and oj^ — 3 X 10~^, in Planck units. We will use these 
values below when we perform numerical integrations. 



with Wq given by 



bibi 



Wo = /Af/o-27^A^5:^(l 



(42) 



The remarkable fact is that all cross terms amongst the 
fields vanish, and Eq. ([40]) is simply the stress energy ten- 
sor for a collection of local (but in general complex) fields 
LP, = with B} = 74Aa[l-4^2(^^_4^^2)jgaa;?^ 

This has been discussed in Q, and a very similar anal- 
ysis, starting with the Lagrangian rather than the field 
equations, was presented in We note that previous 
works which considered multiple roots of the character- 
istic equation took c = in Eq. (|2ip . here we have also 
included the linear term in the potential. 



B. Multiple fields 

In general the characteristic equation ((26|) does not 
have a single root, but as we have seen, a set of roots 
L0I. Moreover, since the non-local field equation ((25|) is 
linear, the most general solution to this equation is given 
by 



(36) 



where (pi are fields satisfying individually 

= -^}i>, - h , (37) 

where we now require the constants to satisfy the ex- 
pression 



E' 



(38) 



ujf m? — 1 
We see that for this general solution we have 

□,^ = -^(c^2^, + 6.) , (39) 

i 

and inserting this into the stress-energy tensor ^ and 
performing the integrations, we arrive at: 



rrilocal 



-^d^ipid^ipi . 



(40) 



C. Quintoms from pairs of complex fields 

For real roots to the characteristic equations we have 
seen that we can rewrite the resulting local fields as 
canonically normalized or phantom fields. However, com- 
plex roots, when viewed in isolation, would lead to com- 
plex local fields and an unphysical evolution (i.e., H and 
consequently the scale factor a would be complex). It 
is important to understand, however, that they can be 
combined in such a way as to produce real fields [S]. 
This follows because the complex roots of Eq. al- 
ways come in pairs which are the complex conjugates of 
one another. For every root ujf, it can readily be verified 
that its complex conjugate (ujf)* is also a root, moreover, 
the associated bi and b* of these roots can also be taken 
to be the complex conjugates of one another. For con- 
venience we relabel pairs of complex roots such that 
is the complex conjugate of o;^^ where m is a positive 
integer, and if real roots exists we label them ui^ (note 
that for certain parameter choices m = n, but this is not 
generically true, for example, in the case ^ there 
can exist a complex root even for n = Oorri = — 1). 

Considering the pair of complex conjugate roots, 
Eq. (j25p implies in turn that ipm and pi^m are also the 
complex conjugates of one another, under the assump- 
tion that the initial conditions from which the fields (pm 
follow are the complex conjugates of the initial condi- 
tions for the fields (p-m- The combination ipm + ip-m is 
therefore real and satisfies local field equations. To un- 
derstand the dynamics further, let us split the complex 
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fields into their real and imaginary parts thus: 

V^m = -^iXm+icrm), (43) 
(P-rn = -^(Xm-icr„), (44) 

where both Xm and <7m are real valued fields. The real 
field which follows from adding the conjugate fields then 
simply becomes (pm + P-m — V^Xm- We also split the 
roots such that, = + Rewriting the local 

stress energy tensor PO)) using this splitting, we find that 
it becomes 



and 



P 



1 

2 ^ 



rrilocal 



r 

\ X! i'^'^Xnid'^Xni - dad.nd" a,n) 

m 

r 777 

(45) 

where the general scalar potential is 

W{ipo, Xm, CTm) = Wo - ^ Sr ^fer<y3r - ^^rV'r^ - 

m 

m 

(46) 

where we have defined bm = [p + i<i)l"/^- We can also 
write 



-Uollfr - br , 



(49) 

are the total energy density and pressure, respectively. 

These field equations can easily be solved for nu- 
merically, so long as we include only a finite number 
of fields in our sum. This is consistent since setting 
Xm = Cm = Xm = o'm — initially, and implies that this 
state will be maintained indefinitely so long as &m = 0, 
hence, we are allowed to make only a finite number of 
quintoms dynamically active. Specifying initial data to 
solve the linear non-local cosmology as an initial value 
problem, therefore, reduces to specifying the initial mag- 
nitude and velocity of each field in our infinite set. 

At this point it is worth recalling a recent analysis in 
Ref. [l^, which sought to understand the initial condi- 
tions required to solve linear non-local theories as an ini- 
tial value problem. Though the analysis presented there 
is more general than the one given here (but less com- 
plete in cosmological settings) , it is clear that the two are 
equivalent as it was found there that each pole of a suit- 
ably defined propagator (equivalent to the zeros of our 
characteristic equation), requires two initial conditions 
(equivalent to our initial field position and velocity). 



V. DIFFUSING LINEAR NON-LOCAL 
COSMOLOGY 

For our linear theory the diffusion-like equation can be 
solved in time once we specify initial functions ^(ti,r) 
and <i'(ti,r) which satisfy the boundary condition (|23p . 
But which class of functions are acceptable for the initial 
data? We found in the previous section that when only 
one field is present, the general form of <I>(t, r), Eq. (|24p . 
is 



□o-m = -/3mXm - al-^CTm - <7m ■ (47) 

Considering G^^ = SttT]^^''^ / together with Eqs. gT]) 
we see that we have succeeded in writing the most general 
local field equations in terms of purely real fields. The 
field X is a regular canonically normalized field whereas a 
is a phantom as its kinetic energy comes with the wrong 
sign. The most general set of local fields which solve the 
non-local equations are therefore given by at most two 
real canonically normalized fields or phantom fields, and 
an infinite set of quintoms. 

So far our work on the correspondence construction has 
been completely general, but as our primary interest is 
in cosmological dynamics, we once again specialize to the 
FRW metric. The field equations become the Friedmann 
equation pH)) and the Raychaudhuri equation where 



^ ( 



1 



(50) 



Hence, a suitable initial profile can simply be taken to be 



$(ii,r) = k-Ce° 



(51) 



W{ipr,Xr,i,(7m) 



(48) 



where k ~ —bjuP'^ C is an arbitrary constant and attP' is 
a solution of the characteristic equation (|26|) . We have 
already studied the possible values for and we found 
that they can have at most two real values and an infinite 
number of complex values. We have also learnt that for 
each complex root, the complex conjugate value is also an 
acceptable root, hence, we can add or subtract C e"" 
and (Ce"" '')* in order that ^(t\,r) and <i>(ii,r) are real. 
Therefore, the initial data can be a linear combination of 
the following: 

$i(ti,r) = 2fcj^-2|C'|e""'«''cos(6'-|-au;?r) , 
$2(ti,r) = 2fc/-2|C|e"'^«''sin(6l + aw|r) . (52) 
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where the indices R and / represent the real and imagi- 
nary parts of (jj^ and k. To confirm that these functions 
do indeed represent suitable initial data, we must calcu- 
late all the infinite constraints implied by the boundary 
condition ((23)) . Let us do this for the linearised model 
defined by ^ o, V{i;) = --iP^/A and V = 1 + (c.f. 
Eg. p^ l. Because of linearity this is simple and we find 



where ^{ti,r) is one of the functions given in Eq. (|52p . 
One can verify that this trial function approximately sat- 
isfies the boundary conditions (fT9|) and that the closer to 
the expansion point A, the better the approximation, and 
therefore it will suit as an initial condition for a numerical 
evolution. 



$(<i,l) = 3$(ti,0)^ 



9$(ti,r) 



Or 



9$(ti,r) 



dr 



r=0 



r=0 



(53) 



and similar for higher order derivatives. It is straightfor- 
ward to verify that the trial functions ((52)) satisfy these 
boundary conditions. One can easily verify that even for 
the more complex cases when ^ 0, the corresponding 
boundary conditions are still satisfied to all orders in the 
linear approximation. 

Before we move to the non-linear non-local case it is 
worth pointing out that for constant H the diffusion 
equation with linear boundary conditions can be solved 
exactly. Here we give the most general solution that sat- 
isfies the boundary condition ((23)) for the case of 6 = 0, 



E 



A. 



[cos(aa;|r) cos(/ii) 



— sm{aiJjr) sin(/it)] 

+ B^^e^-'+°"^> [cos{aLo]r) cos{l2t) 

- sin(au;f r) sin(/2t)] } , (54) 

where A^i and B^^2 are arbitrary constants and 
R± ^ (l/2)Re(-3iJ± V9i?2T^) and I± ^ 
(l/2)Im (-3iJ ± \/9i?2 + aa;2) . This highlights the ill- 
posed nature of our problem, since it can be verified that 
as Loj increases so that the solution becomes more os- 
cillatory in r (which generically occurs for roots to the 
characteristic equation with larger and larger n values) 
the growth rate of the growing mode also increases. 



VI. DIFFUSING NON-LINEAR NON-LOCAL 
COSMOLOGY 

Like in the linear case, if we are to solve the full non- 
linear equations as an initial value problem using the 
diffusion-like equation, then we must specify an initial 
function which satisfies all the infinite conditions that 
we started to write down for the = 0, V{ili) = —if;^ /A 
case in Eq. ((TH)) . In general this is very difficult, how- 
ever, we know that when close to the point about which 
we made our linearised construction, an approximate ini- 
tial function takes the form 



*(ti,r) = A + $(ti,r) 



(55) 



VII. APPLICATIONS AND NUMERICAL 
INTEGRATION 

We have implemented the procedure described in the 
previous sections to compute the time evolution of a field 
i> for two string models. They have in common that the 
general form of the potential in action ([3]) can be written 
as 



V{^) = A 



1 



p+l 



p+i 



(56) 



Hence, the parameters c and m? in the linearised poten- 
tial (dll) read 



= AP -A. 



^pAP- 



(57) 



A. p-adic string 



In the p-adic string A = and the parameters in action 
([3]) correspond to A = 1, = 0, and 



Inp 

2m? 
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1 



(58) 



where gs and are the string coupling and string scale, 
respectively. Expanding the action ([3]) to first order in 
□ we can see that it is equivalent to the action of a field 
with a non-canonical kinetic term 



S4, — 



(59) 



where the effective potential is Vcs{i/j) — 7'*(V^(V')+V'^/2) 
which we illustrate in Fig. (2) for p = 3. Including the 
higher order terms in the derivatives does not change 
the potential part of the action, but affects how the field 
evolves on this potential via the additional kinetic terms. 

In our numerical integration of the equations of mo- 
tion, we have set the initial data by expanding the so- 
lution near the hilltop, A = 1 — e, and have therefore 
fixed 



^-(ti,r) l-ee"'^H''cos(au;?r) 



(60) 
(61) 



= 10- 



is a root given by Eq. 

-12^ 



with e 

ln(j3(l - e)P"^)/(9 X 10-12) ^ith p = 3, and 7^ = 1.48 x 
10~^ in Planck units. We are now ready to evolve the 
diffusion equation ()15p . We have implemented a number 
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FIG. 2: The form of the effective potential for tp in the p-adic 
string with p = 3. The sohd hne and dashed hnes represent 
the full potential and the quadratically expanded potential 
around the hilltop, respectively. 



of finite difference scfiemes and employed tlie method of 
lines, in which the r direction is discretised and ordinary 
differential equation routines were used to integrate the 
resulting equations in t. The results are robust to the 
methods chosen. 

As we mentioned before in Sec. If 11 C[ the numerical 
integration of the diffusion equation is plagued with in- 
stabilities. This is particularly true when H is small or 
negative, since we saw that large H decreases the rate 
at which high frequency modes grow. We therefore re- 
quire to smooth the profile ^{t, r) in r at each time step 
as it evolves, in order to remove the spurious oscillations 
resulting from numerical errors. These usually enter the 
evolution on the scale of the finite grid space in the r 
direction as this represents the highest frequencies per- 
mitted. We perform this smoothing in a number of ways 
in order to ensure the results are not overly dependent 
upon it. Since oscillations appear at high frequencies, one 
possibility is to take an average over a number of points 
around each point in r, in order to smooth the function 
on small scales but leave it unchanged on large scales. 
This can prove difficult to do fairly at points close to the 
boundaries so a procedure with similar results is to take 
a subset of points r*, evaluate ^'(t,r*) at those points 
and perform a spline fitting to reconstruct the data in- 
between. This procedure proves to work very well until 
H becomes negligible. Another possibility, given what we 
have learnt about the form of the profile and the proper- 
ties of the roots to the characteristic equation, is to fit the 
profile to a combination of Eqs. (|52p . Moreover, it can be 
verified that the real parts of the roots vary substantially 
between = 1 and ■)/' = 0, however, the imaginary parts 
are a constant in the case of = 0. Considering this, 
we find that the profile at a given time step can be well 



fitted by 

*fit(t, r) = ci + C2 cl cos(c4 -I- atujr) , (62) 

where ci, C2, C3 and C4 are constants to be determined 
at each time step. Again this theoretically motivated 
method proves highly successful. We have also tried fit- 
ting at each time step to other functions with smooth 
properties such as polynomials. Though with less suc- 
cess, the results show significant improvement over no 
smoothing at all. All these method are forms of regular- 
ization, which generically one must do when dealing with 
ill-posed problems. In fact, our methods are the most 
naive and simplest possibilities, though they appear to 
work well in the setting presented in this article. In fu- 
ture work we hope to return to the issue of regularization 
methods and utilize more sophisticated techniques. 

We use the two methods which work best (the spline 
and the theoretically motivated fitting) to cross check the 
results and we indeed find that they provide the same 
solutions in their respective region of validity. More- 
over we keep track of the accumulated difference S = 
numerical ^'fitl and we choose the size of the 
grid in r and the time step such that this quantity is 
small. As a final check on the validity of our solutions, 
we substitute the solution we find back into a finite dif- 
ference scheme (which can differ from the one we use to 
generate results), and check the finite difference equation 
is satisfied for all r and t. 

The results of the numerical integrations are shown in 
Fig. [3] where the time variable was rescaled with respect 
to the true cosmic time hy t —^ t/^/a for ease of com- 
putation. Setting the initial profile using the n = root 
to the characteristic equation (i.e. Wj = 0), we see that 
the evolution resulting from the linear expansion of the 
potential around the hilltop leads the Hubble rate to van- 
ish when, in the language of the canonically normalized 
field, the potential W{(p) cancels out the contribution 
from the kinetic term (p'^ /2. Once H becomes zero, a 
recoUapse occurs and we do not follow the dynamics any 
further. Evolving the field in the full non-linear poten- 
tial shows that the field settles at the minimum at V' = 
performing damped oscillations and the Hubble rate, H, 
slowly decreases with time. 

This behavior is to be expected because if we make a 
local approximation to the dynamics by expanding about 
the hilltop t/j = 1 the result for the n = root to the char- 
acteristic equation is that the dynamics are equivalent to 
those of a local canonically normalized field evolving on a 
hill-top potential (this is the result of 'g']). On the other 
hand, if we make an expansion at a small value of tp (note 
that though no expansion is possible at '0 = 0, we can 
make the expansion arbitrarily close to "0 ~ 0); the dy- 
namics for the n = root are equivalent to a canonically 
normalized field evolving in a potential minimum made 
up of a quadratic term together with a small linear piece. 
The expected dynamics is for the field to oscillate about 
the minimum of this potential and to decay towards it. 
The non-linear numerical integration appears to interpo- 
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FIG. 3; Upper left panel: Evolution of the field ip in the p- 
adic string, = 0, p = 3, where ip is & single real field, n = 0, 
and initially A = 1 — t. The solid and dashed lines represent 
the non-linear and linear evolution, respectively. Upper right 
panel: The corresponding evolution of the Hubble rate. Lower 
left panel: The number of e-folds before the end of inflation. 




t 



FIG. 4: Upper left panel: Evolution of the field %p in the p- 
adic string, = 0, p = 3, where t/j is a single real field, n = 0, 
and initially A= 1 + e. The solid and dashed lines represent 
the non-linear and linear evolution, respectively. Upper right 
panel: The corresponding evolution of the Hubble rate. Lower 
left panel: The number of e-folds before the end of infiation. 



late between these two behaviors. In Fig. [31 we also show 
the number of e-folds before the end of inflation, N , and 
we can clearly see that when A'^ > 60, the linear and non- 
linear evolutions are practically indistinguishable. This 
observation suggests that the computation of the scalar 
power spectrum using only the linear approximation in 
Ref. ^] is sufficiently accurate. It is also interesting to 
compare our solution with the solution recently found 
for the light-like p-adic tachyon in a linear dilaton back- 
ground which also evolves to the minimum of its potential 

In Fig. 3] we show the evolution of the same quanti- 
ties as in Fig. [3] but with the field starting its evolution 
ai A — 1 -f e. Clearly the linear evolution is identical 
to the previous situation but the non-linear dynamics is 
different as there is no minimum to stabilize the field on 
the right hand side of the hilltop. Though the solution is 
not of particular interest, it highlights the utility of our 
method as opposed to a boundary value formulation as in 
that approach the field is only allowed to evolve between 
two extrema. 

With n — \, the first quintom, the field performs grow- 
ing oscillations around the hilltop as can be seen in Fig. [5l 
These oscillations are mimicked by the Hubble rate and 
eventually they are so large that cause this quantity to 
vanish. This occurs very near the hill top, hence, the 
linear and non-linear evolutions are practically indistin- 
guishable. Once H goes to zero we do not follow the 
evolution any further as the universe recollapses at this 
point. Although a re-expansion is a possibility, such 
cyclic dynamics is more difficult to deal with as negative 
values of H lead to even greater numerically instabili- 
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FIG. 5: Upper left panel: Evolution of the field tp in the 
p-adic string, = 0, p = 3, where Tp is the first quintom, 
n = 1. The solid and dashed lines represent the non-linear 
and linear evolution, respectively. Upper right panel: The 
corresponding evolution of the Hubble rate. Lower left panel: 
The number of e-folds before the end of infiation. 



ties, and in any case, they are beyond the scope of our 
investigation. 

We would like to emphasize once more that though 
we used a small value of p in our examples, the effect of 
higher order derivatives in the equation of motion is still 
present and does play an important role. Let us make 
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this more clear. If we had started from action (|59|) that 
only considers the first order term in □, and expanding 
using xp — 1 + (f), we would have obtained the equation of 
motion 



2 

— ( 

a 



(63) 



However, the local equation of motion (|25)) that follows 
from the full linearised non-local equation of motion (PO]) 
gives instead 



□(/)=- 



(64) 



which, though similar for small p, is clearly different from 
Eq. ([63]) and leads to great disparity in the results for 
large p. 



B. Cubic superstring field theory (CSSFT) 

In the CSSFT case we have A = -1, = 0.95, p = 3, 
and 



1 



a 



1 



(65) 



4mf 

Expanding the action ^ to first order in □ , we obtain 

(66) 

Hence, for ^■^ > 1/4, similarly to what happens in the 
p-adic case, we also have a non-phantom field with non- 
trivial kinetic terms. The effective potential is in this 
case Veff (V') = —l'^{V{i^) + V'^/2), where A is a constant 
related to the D-brane tension and we set its value such 
that V^ff(l) = 0. This potential is shown in Fig.[6l 

The numerical procedure to solve the diffusion equa- 
tion in this model follows closely the method described 
for the p-adic string, but since now the hilltop is at ^/i = 
we set A = e and hence, we specify the initial data as 



■^{tur) = ee""«''cos(at^|r) 
^'(ti,r) = 0, 



(67) 
(68) 



with e = 10 is a root given by Eq. ([^ . a = 1/(36 x 

10^12 ^4 ^ 48 ^ ;^Q-9 Planck units. For ^ 

0, however, the imaginary parts of do vary with the 
expansion point A but only by a significant amount when 
n = Oorn = — 1. In these cases, the imaginary parts only 
exist for |^| > 0.7 and correspond to small frequencies 
as auj] < 1, hence, one period is not completed between 
< r < 1. For n > 1 and n < —1, the frequency of 
the profile is nearly independent of A. In this model, 
therefore, we also use as a fitting smoothing function the 
profile 

The results of integrating the diffusion equation nu- 
merically are shown in Fig. [71 Starting the field near 




FIG. 6: The form of the effective potential for tp in the CSSFT 
case. The solid line and dashed lines represent the full po- 
tential and the quadratically expanded potential around the 
hilltop, respectively. 



Or 

-0.5 

-1 

80 r 

60 
^40 
20 



1.5 

1 ' 
0.5 



X 10" 



H 



50 100 150 




50 100 150 




50 100 150 



t 



FIG. 7: Upper left panel: Evolution of the field tp in the 
CSSFT case, = 0.95, p = 3 where tp is originally a sin- 
gle real field, n — 0, but develops a quintom for \tp\ > 0.7. 
The solid and dashed lines represent the non-linear and linear 
evolution, respectively. Upper right panel: The corresponding 
evolution of the Hubble rate. Lower left panel: The number 
of e-folds before the end of inflation. 



the hilltop corresponds, as we have seen, to a real field. 
Once it has been displaced beyond \ip\ w 0.7 it develops a 
quintom as auj'^ becomes imaginary. We indeed sec that 
the Hubble rate performs an oscillation when the field is 
near |^| = 1 and eventually collapses to H = 0. As in the 
p-adic case, the non-linear analysis interpolates between 
the linear expansion behavior about the hilltop and the 
linear expansion evolution about the minimum, where in 
this case, only quintoms are found. 
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FIG. 8: Upper panel: Evolution of the field tp in the CSSFT 
case, = 0.95, p = 3 where ip is a quintom with n = 1. 
The solid and dashed lines represent the linear and non-linear 
evolution, respectively. Lower panel: The corresponding evo- 
lution of the Hubble rate. 



As we have seen in Section IIVI when = exactly, 
there are no quintoms allowed as the only root is real. 
But slightly away from the hilltop, however, there can 
exist quintoms as the root can be complex (see Eq. (|29p ). 
The evolution for the quintom with n — 1, is shown in 
Fig. [8l Its behavior is somewhat surprising given that 
it is fairly distinct from the evolution of the quintoms in 
the p-adic string (in Fig.O- Instead of oscillating around 
the hilltop with growing amplitude, the amplitude of the 
field decays around the hilltop. This can be explained by 
noticing that the real part of is negative and becomes 
larger in absolute value when the field is near the hilltop 
whereas its imaginary part is small and nearly constant. 
This implies that both modes in solution ([51)1 are decay- 
ing. In this example, the field gets trapped at V' = 
and consequently the expansion of the universe becomes 
closer to de Sitter. Though this quintom is a viable way 
of obtaining accelerated expansion, an additional mech- 
anism to exit inflation and enter the radiation epoch is 
required to complete the scenario. 



VIII. AN ASIDE ON PERTURBATIONS 

Our analysis of non-local dynamics can readily be ex- 
tended to perturbations about the homogeneous back- 
ground evolution. We leave a full analysis to future work, 
and simply make some observations here. 

First, we note that for linear non-local theories the 
localization technique works equally well when a pertur- 
bation is added to the background field and to the metric. 
This is valid because the technique uses an unspecified 
□ operator. Then the perturbed linear non-local case 
is simply equivalent to a perturbed infinite field local 



cosmology. The diffusion-like equation approach to the 
perturbed linear non-local theory must therefore also be 
tractable in this case since it is equivalent to the local- 
ization technique. For the non-linear case the situation 
is more involved but a similar idea to that employed for 
the background is possible. That is, we can use the ini- 
tial profiles determined in the linear case as approximate 
initial conditions for the non-linear case. 

As a first approximation one might consider perturb- 
ing the non-local scalar field and neglecting the inevitable 
back-reaction on the metric keeping it unperturbed. Sep- 
arating out the unperturbed equations, and decomposing 
the field perturbation into Fourier modes, the non-local 
scalar field equation for a particular mode takes the form 



Sipk 



(69) 



where k labels the wave number of the mode, and where 
□fc is 

fc2 



D - 3ff^ 



(70) 



H and a are simply the Hubble rate and scale factor of 
the unperturbed universe respectively. 

The first thing to note is that this equation is linear, 
but the situation is not as straightforward as in the linear 
unperturbed case as the effective mass (d'^V/dilj^) is now 
time dependent. This means that the correspondence 
technique does not work. We could approximate the ef- 
fective mass as a constant for a short period of evolution, 
for example, as a mode exits the cosmological horizon but 
this procedure is bond to be inaccurate. An alternative 
is to use the diffusion-like equation technique. We can 
introduce the new field 



(5*(t, r) 
which satisfies the equation 



SMt): 



1 9(5* 



(71) 



(72) 



a dr 

For the potential V{ip) = -ip^/A and = o, the field 
equation for the perturbation translates into the bound- 
ary condition 

(5*(t, 1) = 3*(t,0)2,5*(t,0), (73) 

which is linear but time dependent. As for the unper- 
turbed field, we must specify ^'(ti,r) and 5'(ti,r) which 
must also satisfy the boundary conditions (|19p . if we are 
to solve the diffusion equation as a initial value problem. 
However, a very similar problem arises as in the non- 
linear background case. That is, our single boundary 
condition in the r direction actually imposes an infinite 
set of boundary condition. Indeed, the second boundary 
condition with 5'(<i,r) = (5'I'(ii,r) ~ reads 



S^*(ii,r) 
dr 



r=l 



6*(<i,0) 
+3*(ii,r)2 



9*(ti, 



a(5*(ii,r) 



dr 



'5*(ti,0) 
.(74) 
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We can only satisfy these approximately as we did for the 
non-linear unperturbed case, and then use the diffusion- 
like equation to evolve forward in time. While we leave 
such an analysis to future work, we mention it here to 
highlight the utility of the initial value formulation of 
the diffusion equation method when compared with the 
boundary value formulation. Of course one could also 
obtain a solution using the boundary value formulation, 
but for that procedure 6'^ must tend to a constant at 
±00 fixed by Eq. For example, for the p-adic case 

we must have Sif'k = S-i/j^iJ^fe at t = ±00 as all the 
derivatives vanish there. Since ip can only take values 
■0 = ±1,0 at t = ±00, it is clear that Sipk = 0. Hence 
the boundary value formulation can only generate solu- 
tions in which the perturbation decays asymptotically. 
However, in a study of perturbations produced during 
non-local inflation for example, we would like to put a 
particular /c-modc of the perturbation into a state simi- 
lar to the Bunch-Davis vacuum initially and evolve it for- 
ward in time until it exits the horizon. Repeating this for 
many modes we can determine the spectrum of pertur- 
bations produced. This utility is a major motivation for 
studying the diffusion-like pde technique from the point 
of view of an initial value problem. 

IX. CONCLUSIONS AND FUTURE 
DIRECTIONS 

In this work, our major concern has been whether non- 
local cosmologies can be solved as an initial value prob- 
lem by solving a related partial differential equation as an 
initial value problem. We have found that this is indeed 
possible but that there are technical difficulties, both in 
understanding how to specify suitable initial data and in 
overcoming the ill-posedness of the pde. 

The first difficultly is easily resolved in linear models 
where separation of variables leads to an interpretation 
of the non-local cosmology in terms of a local one with 
an infinite number of fields. We showed, intriguingly, 
that most of these fields are in the form of quintoms (one 
standard field and a coupled phantom field). Moreover, 
the separation of variables technique led to a simple un- 
derstanding of suitable initial data for the diffusion-like 
equation. In non-linear cases, we showed that linearis- 
ing the equations allowed approximate initial data to be 
specified for a subsequent non-linear evolution. Extend- 
ing previous analysis, we showed that this linearisation 
can be preformed about regions of the potential other 
than a maximum or minimum. 

The second difficultly we overcame by using a sim- 



ple regularization technique, where at each time step we 
smoothed the profile of the solution in the r direction 
to eliminate spurious oscillations which would otherwise 
swamp the solution. The smoothing which worked best 
took advantage of the theoretical knowledge about the 
probable form of the profile (from the linearised anal- 
ysis), though we stress that the free parameters in the 
profile were fitted from the data rather than imposed. 

Previous work on inflation in non-local theories has 
been carried out by a number of authors. A major conclu- 
sion of our work is that the method employed in Ref. Q 
(i.e. the linearisation of the equations about the hill-top, 
considering the evolution of only the real canonically nor- 
malized local field) works very well for the background 
evolution when observationally relevant scales are exiting 
the cosmological horizon in the cases we have studied. 
Though we note that the evolution can be very differ- 
ent if other local fields (the quintoms) are permitted to 
be present initially. Such a localization method allows 
known results for a canonically normalized field to be 
employed, and hence the power spectrum can be calcu- 
lated for the non-local theory. Of course, before making 
a firm statement we would have to perform the anal- 
ysis presented in this article for perturbations as well 
as for the background evolution. We sketched how this 
could be achieved in Section IVIIII and we will return to 
this issue in future work. Many questions should be rig- 
orously addressed, for example: In a non-local theory, 
is the comoving curvature perturbation still conserved 
on super-horizon scales even when the evolution is non- 
linear? This needs to be known before the previous re- 
sults (and in particular claims of a large non-gaussianity 
in p-adic inflation when p is large 7]) can be confirmed 
or falsified as they all relay on it. 

Since several authors have also considered non-local 
tachyons as the source of dark energy in the universe 

, another obvious application of our methodology is to 
include a background fluid into these models and see how 
the field behaves during radiation and mater domination. 
Only when this is known can a firm statement on whether 
a non-local field is a dark energy candidate can be made. 
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